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Abstract 

Coexistence of antiferromagnetism and d-wave superconductivity within a Cu02 plane was re- 
cently observed in a wide doping region for multilayer high-temperature cuprate superconductors. 
We find that the experimental phase diagram is well reproduced in the slave-boson mean-field 
scheme of the two-dimensional t-J model by including antiferromagnetic order. We argue that 
weak three dimensionality coming from a multilayer structure is sufficient to stabilize antiferro- 
magnetic order and its coexistence with superconductivity. 

PACS numbers: 74.20.Mn, 74.72.-h, 74.25.Ha, 71.10.Fd 



I. INTRODUCTION 



High-temperature superconductivity (SC) is one of the most fascinating topics in physics. 
The discovery of the magnesium diboride in 2001 (Ref. [l|) immediately generated very in- 
tensive studies. The recent discovery of ferropnictides^ also attracts tremendous attention, 
opening a new branch in condensed matter physics. It is however still only cuprate supercon- 
ductors especially multilayer cuprates such as HgBa2Ca„„iCunOy and TlBa2Ca„„iCu„0j^ 
with n > 3 that achieve superconducting transition temperature (Tc) more than 100 K 
(Ref. 13|); we mean by "multilayer" three or more layers in a unit cell in this paper. For such 
multilayer cuprates, the previous theoretical studies^"- payed attention to the description of 
a superconducting phase coherence between the layers, including a possible charge imbal- 
ance between the layers in a unit cell. In those studies there was a tacit assumption that the 
property of each Cu02 plane is essentially the same as that in single- and bi-layer cuprates 
in the sense that AF is realized for very low doping and is replaced by SC for moderate 
doping without the coexistence with AF within a Cu02 plane. 

However, the recent NMR measurements for multilayer cuprate superconductors^i"- re- 
vealed the phenomena very different from those in single- and bi-layer cuprates: antiferro- 
magnetism (AF) in the Mott insulator survives up to rather high carrier doping and coexists 
with a superconducting state. The coexistence was due to not a proximity effect between 
the layers, but a phase transition within a Cu02 plane at low temperatures. 

This observation sharply contrasts with a widely accepted viewpoint that AF in the 
Mott insulator is rapidly suppressed by a tiny amount of carrier doping and the cZ-wave 
superconducting state is instead stabihzed in a wide doping region, typically in 5 — 30%. 
The latter viewpoint is based on the data for single- and bi-layer cuprates such as La- and 
Y-based compounds, for which a huge number of studies have been performed last about 25 
years in a systematic way. Compared with those materials, multilayer cuprates are much 
less investigated. They contain completely flat Cu02 planes with a perfect square lattice 
and are known to be free from disorder, in contrast to La- and Y-based cuprates. They 
also achieve higher than that for single- and bi-layer cuprates.^ In this sense multilayer 
cuprates can be viewed as an ideal system to study the mechanism of high Tc. Nonetheless, 
the basic theoretical framework for multilayer cuprates has not been identified so far. 

Moreover the origin of the pseudogap^^ is still a central issue on high-temperature SC 
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and a clue to resolve it is highly desired. Therefore it is of great importance to shed light 
on an issue of the pseudogap in multilayer cuprates from a theoretical point of view, which 
may in turn provide a crucial insight into the pseudogap in single- and bi-layer cuprates. 

In this paper, we explore the basic theoretical framework which captures the essential 
features recently reported for the multilayer cuprates." We invoke Anderson's resonating- 
valence-bond (RVB) scenario,— that is, the undoped Cu02 plane is assumed to be in the RVB 
spin liquid state and the preexisting spin singlet pairs can become charged superconducting 
pairs once they are mobile by carrier doping. This idea was well described in terms of 
the two-dimensional (2D) t-J model. In particular, the slave-boson mean-field theory^^^^^ 
and the gauge theory,^^ which takes low-energy fluctuations around the mean-fields into 
account, turned out to capture many important properties of single- and bi-layer cuprate 
superconductors.—"— 

In the standard RVB framework, AF is assumed to be strongly fluctuating, not to be 
ordered. However, antiferromagnetic order can be easily stabilized in the presence of weak 
three dimensionality coming from a multilayer structure. This effect can be incorporated 
phenomenologically by including AF as a possible mean-field order parameter within a single- 
layer model. Such a calculation was performed in the slave-boson mean-field scheme^*^ and 
variational Monte Carlo^i^i for the 2D t-J model in a different context, showing that AF 
extends to a high doping rate and coexists with SC. However, the coexistence obtained 
previously was found to be substantially suppressed by the presence of a long-range hopping 
amplitude, leading to a pure antiferromagnetic state.— Moreover, the previous theoretical 
results^*^ typically exhibited reentrant behavior of the critical temperature of AF at low 
temperatures, which was not observed experimentally.- 

We perform the slave-boson mean-field analysis of the t-J model by including antifer- 
romagnetic order. While we analyze the 2D model and do not take multilayer degrees of 
freedom into account, we consider that a special feature of multilayer cuprates is included 
phenomenologically by allowing AF order as a possible mean field in our analysis. In this 
framework, important features of the phase diagram for multilayer cuprates are obtained.— 
We argue that the essential difference between multilayer cuprates and single- and bi-layer 
cuprates lies in the presence of weak three dimensionality and thus the phase diagram of the 
former is well described by both RVB and AF, while that of the latter is simply by the RVB 
as already investigated. In contrast to single- and bi-layer cuprates, the pseudogap may be 
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substantially diminished in the heavily underdoped region of multilayer cuprates. 

II. MODEL AND FORMALISM 

We employ the 2D t-J model on a square lattice 

"H = - ^ tijc\„Cja + J ^ Sj ■ , (1) 

where the transfer integrals are finite for the first- (t), second- (t'), and third-nearest 
neighbor bonds (t"), and vanish otherwise. J(> 0) is the antiferromagnetic superexchange 
interaction and («, j) denotes the nearest neighbor bonds. is the electron operator in 
the Fock space without double occupancy and we treat this condition using the slave-boson 
method^^ by writing Ci^ = h\fia under the local constraint fiafi<T + = 1 at every i 
site. Here /jo- {hi) is a fermion (boson) operator that carries spin a (charge e); the fermions 
(bosons) are frequently referred to as spinous (holons). The spin operator is expressed as 
^i = \ Hap fia^a^fip ^ith CT bciug the Pauli matrices. 

Hamiltonian ([T]) is decoupled by introducing the following order parameters:—!^ the stag- 
gered magnetization m = \{f-^fi^ — /^/ii)^**^''"' with Q = (vr, vr); the bond order parameter 
for spinous and holons, {J2a fLfj^)^ '^^ denote x = {J2a flafjo) ^^e nearest neigh- 

bor bond; the singlet RVB paring A,- = {fi^fi+ri — fiifi+rt) with t = x,y. Here we assume 
that all these expectation values are real and independent of i. We can show that the dx2-y2- 
wave pairing state is the most stable, i.e., A^, = —Ay = Aq. Although the bosons are not 
condensed in the present mean- field scheme at finite temperature (T), they are almost con- 
densed at low T and for finite carrier doping 5(> 0.02).— Hence we approximate (6) ~ ^/S 
and (blbj) ^ 6. In principle, the so-called vr-triplet pairing state can emerge for a state with 
m 7^ and Aq 7^ 0,^^ but turns out not to be stable in our model. Hence the free energy 
per lattice site is computed as 

+ ^{x' + Al)+2Jm'-l^5, (2) 

where = a/ (^k )^ + spinon's band dispersion in the presence of m and Ag; 

= C ± D^, = V(C)' + (2^^)^ ^k = (^k ± ^k+Q)/2, ^k = -2{t6 + |Jx)(cos A;,. + 
cos ky) — At' 6 cos kx cos ky — 2t" 6 {cos 2kx + cos 2ky), and Ak = — | JAo(cos /c^: — cos fcy); /i and 
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denote the chemical potential and the total number of lattice sites, respectively; the sum 
of momentum is taken over the magnetic Brillouin zone I hx \ ~t~ | ky | ^ tt . 

Because we relax the local constraint to a global one fiafi^) = 1 — ^ and = S 

in the present mean-field theory, our approximation may be reliable as long as electrons are 
in coherent motion. In the present case, since x tends to saturate for 6 > 0.05 [Fig. [Hb)], 
we expect that our approximation is sufficiently reliable in such a region, where most of 
experimental data have been obtained so far. 

To examine a possibility of the incommensurate antiferromagnetic instability, we also 
compute the longitudinal magnetic susceptibility x(q) in the random phase approximation: 
X(q)"^ = Xo(q)"^ + 2J(cosg^ + cosgj^) and 



AN 

k 



tanh ^ - tanh %a tanh ^ + tanh %a 

k,k+q ^ ^ ^ '-'k,k+q 7-1 , 7-1 



(3) 



± _ W 1 _U '^k^k+q + AkAk+q 



Here = a/^^ + A^ is the spinon's band dispersion for m = and the sum of k is taken 
over the region \kx\, \ ky\ < n. 

The material dependence of cuprate superconductors can be taken into account mainly by 
different choices of t' and t".—'— Hence it is naturally expected that we could invoke specific 
values of t' and t" appropriate for multilayer cuprates. This could be achieved for a realistic 
multilayer model, which also contains other parameters such as interlayer hopping integrals, 
interlayer exchange interactions, and site potential yielding a charge imbalance between 
the layers in a unit cell. Given that there is much ambiguity about those parameters and 
that a special feature of multilayer cuprates is considered phenomenologically in the present 
analysis by invoking AF order as a possible mean field, we consider the values of t' and t" 
simply as phenomenological parameters to reproduce various types of the phase diagram. 



III. RESULTS 



We determine the mean fields by minimizing the free energy Eq. ([21) and obtain the phase 
diagram in the plane of 6 and T. The result for t/J = 4, t' /t = 0.12, and t"/t = —0.06, 
for which an electron-like Fermi surface is realized in a normal state [Fig. lU^c)], is shown 
in Fig. [D^a). The temperature T/v denotes the onset of AF, and (Trvb and T^vb^) 
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the onset of singlet paring in the presence (absence) of AF. At 5 = 0, AF is reahzed and 
no singlet pairing coexists. With carrier doping, the critical temperature of AF gradually 
decreases and survives in a wide doping region up to 6^ ~ 0.17; the 6]\f is the critical doping 
rate of AF at T = 0. AF suppresses the formation of singlet pairing (T^-q < T^yg^), 
but the coexistence is realized at low temperatures. Singlet pairing extends over a wider 
doping region than AF, and a pure (i-wave singlet pairing state is realized at high doping 
rates. The right-hand panel in Fig. [T]^a) magnifies the region around the tetracritical point 
(5tct, ^tct) ~ (0.165, 0.031J), where the four states, AF, SC, their coexistence, and the normal 
states, become identical. The order parameters m, Aq, and x ^i'g shown in Fig. [Il^b) as a 
function of 6 at low T. Although Aq still increases with decreasing 5 in a certain region in 
^ ^ ^AT, we see that Aq is typically suppressed by the presence of AF compared with Aq for 
m = 0. Similarly the magnitude of x is suppressed by the presence of m and becomes zero 
at (5 = 0. The system becomes an antiferromagnetic insulator at 5 = 0. 

A crucial feature of Fig. [T] (a) (see the top right-hand panel) is that in spite of competition 
of AF and singlet pairing, AF extends to a higher doping region beyond the tetracritical 
point in the singlet pairing state. This is the crucial difference from the previous work— 
where T/v exhibits reentrant behavior at low T, leading to 6^ < 5tct- The difference comes 
from our careful choice of the band parameters under the constraint t" = —t'/2 (Ref. |29| ) 
such that the shape of the Fermi surface around 6 = 6tet fulfills the nesting condition of 
q = Q close to the nodal region of the (i-wave pairing gap. Hence, the resulting static spin 
susceptibility x(q) shows a maximum at q = Q and is not suppressed by the onset of singlet 
pairing, leading to 6n > 5tet- 

The careful choice of the band parameters is not a unique way to obtain a phase diagram 
similar to Fig. [Tt^a) when we allow other magnetic orders with q 7^ Q. In Fig. Et^a), we 
show a representative result. We first minimize the free energy Eq. ([2]) with respect to 
the mean fields and Aq. We then obtain the lines of T/v, T^yb^ ^rvb- The 

Tat exhibits reentrant behavior around 6 = 6n ^ 0.15 at low temperatures, in contrast to 
Fig. [U^a). To investigate a possibility of other magnetic orders, we check a wave vector 
at which x(q) diverges. It turns out that a part of the line T/v [solid line with crosses in 
Fig. El^a)] is preempted by the onset of incommensurate antiferromagnetic order T^, where 
x(q) diverges at q 7^ Q. The ordering wave vector q depends strongly on 6 and T. The 
most crucial point in Fig. [2]^a) is that the original reentrant behavior of T/v (solid line with 
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FIG. 1: (Color online) (a) Phase diagram in the plane of 5 and T for t/J = 4, t' /t = 0.12 
and t" /t = —0.06. T/v is the onset temperature of AF, and T^-^ (Trvb and TjlJvB^) ^^s^t of 
singlet pairing in the presence (absence) of AF. The right-hand panel magnifies the region around 
the tetracritical point; the antiferromagnetic and superconducting states are denoted by "m" and 
"Aq", respectively, and their coexistence is by "Ao + m". (b) 5 dependence of the order parameters 
at T = O.OOIJ. (c) Fermi surfaces in the normal state for several choices of 5. 

crosses) is fictitious and instead a true phase boundary is given by T^. We find that these 
results are generic and applicable for various band parameters, which reproduce an electron- 
like Fermi surface that does not cross the magnetic Brillouin-zone boundary + l^yl = 
around 5 = 5tet [see Fig. Et^b)]. In this sense, the phase diagram of Fig. HJa), where AF 
with q = Q is stabilized, should be regarded as a special case requiring a tuning of band 
parameters as already mentioned above. In fact, if we change the value of t/J in Fig. [T]to 
t/J = 3, keeping t'/t and t"/t unchanged, the Fermi surface stays almost the same, but 6tet 
shifts to be a bit larger. Such a small shift is sufficient to degrade the nesting condition 
of q = Q around 6 = 5tct- The resulting T/v exhibits reentrant behavior, which is however 
preempted by T^, similar to Fig. [2]^a). There would also be a possibility that the reentrant 
behavior of T/v shown in Fig. |2]^a) could be preempted by a first order transition to AF with 
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FIG. 2: (Color online) (a) Phase diagram in the plane of S and T for t/J = A, t' = t" = 0. T^^ 
is the onset temperature of incommensurate AF; other notations follow those in Fig. [^a). AF is 
realized in the shaded region, (b) Fermi surfaces. The boundary of the magnetic Brillouin zone is 
denoted by a dotted line. 




FIG. 3: (Color online) (a) Phase diagram in the plane of 5 and T for t/J = 4, t' /t = 0.3, and 
t" /t = —0.15. The notations follow those in Fig. HJ^a). (b) Fermi surfaces. The boundary of the 
magnetic Brillouin zone is denoted by a dotted line. 

q = Q. However we checked that such a possibility does not occur by observing that the 
Landau free energy features a single minimum as a function of m. 

If an electron-like Fermi surface crosses the magnetic Brillouin-zone boundary around 
5 = Stet as shown in Fig. MJo), which may be applicable to electron-doped cuprates,^°'-i T/v 
tends to exhibit a straight line but still features a continuous phase transition [Fig. El^a)]. 
Incommensurate magnetic order is not found to be stabilized. On the other hand, for band 
parameters leading to a hole-like Fermi surface such as that frequently used theoretically for 
Y- and Bi-based cuprates, AF with q = Q is the most stable around 6n and the line of Tn 
exhibits reentrant behavior as already seen in the literature.— 
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IV. DISCUSSION AND CONCLUSION 



Our obtained phase diagrams Figs. [T] and [2] show good agreement with the experimental 
phase diagram.-i^ AF extends over a wide doping region and coexists with SC at low T. 
Not only the value of 6^ but also the magnitude of m is comparable with the experimental 
results.-"- We have invoked an electron-like Fermi surface such as that shown in Figs. [U^c) and 
[2t^b). We expect such a Fermi surface is stabilized as one of Fermi surfaces in multilayer 
cuprates^^ by virtue of strong hybridization between the layers with the dispersion cx 
(cosfc^ — cos ky)"^ (Ref. l29|). The angle-resolved photoemission spectroscopy (ARPES) was 
performed for Bi2Sr2Ca2Cu30io+y— and Ba2Ca3Cu408F2,— and failed to reveal all Fermi 
surfaces. In connection with the NMR experiments,-"- it is desirable to perform the ARPES 
for Hg- and Tl-based cuprates and to test a possible presence of an electron-like Fermi 
surface. It is not clear which phase diagram, Fig. [T](a) or Fig. [2]^a), is more appropriate to 
multilayer cuprates, since the NMR does not directly discriminate different ordering patterns 
of magnetism. Hence it is also desirable to perform neutron scattering measurements for 
the multilayer cuprates to reveal the wave vector of magnetic order. 

The singlet paring formation of spinous is interpreted as the pseudogap in the slave- 
boson formalism in the underdoped region. Since the optimal carrier density 6op of mul- 
tilayer cuprates is situated above the tetracritical point 6tet, ^rvb in ^tet ^ ^ ^op in 
the phase diagram is interpreted as the pseudogap temperature T*. The previous NMR 
measurements^*^ indeed observed the pseudogap behavior of (TiT)~^ in such a doping re- 
gion. Since J is around 100-150 meV, the obtained value of Trvb in Figs. [T] and [2] is small 
compared with the experimental observation.— This discrepancy should be explored by 
including explicitly the multilayer degree of freedom in the present analysis. 

On the other hand, for 6 < 6tct, U(l) gauge fluctuations emerging in the slave-boson 
formalism are expected to be strongly suppressed in the antiferromagnetic state and thus 
spinous and holons tend to confine there.— Therefore we expect that the pseudogap in the 
antiferromagnetic state is substantially diminished for 6 < 6tet and instead the coexistence of 
AF and SC is realized unless the tendency of carrier localization appears at low temperatures. 
Around the tetracritical point it is interesting to clarify both theoretically and experimentally 
how the onset temperature of the pseudogap changes to that of the phase transition to the 
coexistence with decreasing 6. 
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A crucial difference between single-layer and multilayer cuprates lies in the difference of 
antiferromagnetic fluctuations. AF is realized through breaking of continuous symmetry, 
which thus does not occur at a finite T in a pure 2D system.— Hence the presence of 
AF in layered materials is interpreted as coming from weak three dimensionality, which 
is always present in real systems and in general suppresses fluctuations. In single-layer 
cuprates, because of the intrinsic low dimensionality due to a tiny coupling between Cu02 
layers along the c axis, antiferromagnetic fluctuations are expected to be so strong that the 
magnetism is realized only close to the Mott insulator. Moreover some extrinsic effect such as 
randomness may easily hinder long-range antiferromagnetic order. The standard slave-boson 
formalism of the t-J model (without AF) was proposed for such cuprate superconductors.— 
On the other hand, for multilayer cuprates, many layers are already present within a unit 
cell, yielding relatively strong three dimensionality compared with single-layer systems. In 
addition, each Cu02 plane is perfectly flat and free from disorder. These can be the main 
reasons why the present slave-boson mean-field analysis with AF captures essential features 
observed in multilayer cuprates.-"- A natural consequence is that AF would extend over a 
wider doping region by increasing the number of Cu02 planes in a unit cell unless the value 
of t/J varies significantly. This tendency is actually reported in Ref. js]. 

Multilayer cuprate superconductors achieve much higher Tc than single-layer systems. 
Hence the understanding of multilayer cuprates is crucially important to elucidate the mech- 
anism of high-temperature SC. We have argued that weak three dimensionality coming from 
a multilayer structure is sufficient to stabilize AF and that multilayer cuprates can be sys- 
tems described by RVB and AF in the t-J model, suggesting the importance of the local 
antiferromagnetic coupling J. 

Acknowledgments 

We are grateful to H. Mukuda for useful discussions. H.Y. thanks O. K. Andersen, M. 
Fujita, V. Hinkov, A. A. Katanin, W. Metzner, and R. Zeyher for valuable discussions. 



J. Nagamatsu, N. Nakagawa, T. Muranaka, Y. Zenitani, and, J. Akimitsu, Nature 410, 63 
(2001). 



10 



2 Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. Am. Chem. Soc. 130, 3296 (2008). 
^ A. lyo, Y. Tanaka, H. Kito, Y. Kodama, P. M. Shirage, D. D. Shivagan, H. Matsuhata, K. 

Tokiwa, and T. Watanabe, J. Phys. Soc. Jpn. 76, 094711 (2007). 
^ S. Chakravarty, H.-Y. Kce, and K. V51kcr, Nature 428, 53 (2004). 

5 T. A. Zaleski and T. K. Kopec, Phys. Rev. B 71, 014519 (2005). 

6 M. Mori and S. Maekawa, Phys. Rev. Lett. 94, 137003 (2005). 

^ H. Mukuda, Y. Yamaguchi, S. Shimizu, Y. Kitaoka, P. Shirage, and A. lyo, J. Phys. Soc. Jpn. 
77, 124706 (2008). 

^ S. Shimizu, T. Sakaguchi, H. Mukuda, Y. Kitaoka, P. M. Shirage, Y. Kodama, and A. lyo, 

Phys. Rev. B 79, 064505 (2009). 
^ S. Shimizu, H. Mukuda, Y. Kitaoka, H. Kito, Y. Kodama, P. M. Shirage, and A. lyo, J. Phys. 

Soc. Jpn. 78, 064705 (2009). 
1° T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). 

11 P. W. Anderson, Science 235, 1196 (1987). 

12 G. KotUar and J. Liu, Phys. Rev. Lett. 38, 5142 (1988). 

1^ Y. Suzumura, Y. Hasegawa, and H. Fukuyama, J. Phys. Soc. Jpn. 57, 2768 (1988). 
1^ N. Nagaosa and P. A. Lee, Phys. Rev. Lett. 64, 2450 (1990). 

15 p. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 78, 17 (2006). 

16 H. Yamase and W. Metzner, Phys. Rev. B 73, 214517 (2006); H. Yamase, ibid 75, 014514 
(2007). 

1^ M. Ogata and H. Fukuyama, Rep. Prog. Phys. 71, 036501 (2008). 

18 M. Inaba, H. Matsukawa, M. Saitoh, and H. Fukuyama, Physica C 257, 299 (1996). 

19 H. Yamase and H. Kohno, Phys. Rev. B 69, 104526 (2004). 

20 A. Himeda and M. Ogata, Phys. Rev. B 60, R9935 (1999). 

21 C. T. Shih, Y. C. Chen, C. P. Chou, and T. K. Lee, Phys. Rev. B 70, 220502(R) (2004). 

22 A part of the present work was pubhshed as conference proceedings, K. Kuboki, M. Yoneya, 
and H. Yamase, Physica C 470, S163 (2010). 

23 Z. Zou and P. W. Andersen, Phys. Rev. B 37, 627 (1988). 

24 G. C. Psaltakis and E. W. Fenton, J. Phys. C 16, 3913 (1983). 

25 T. Tanamoto, H. Kohno, and H. Fukuyama, J. Phys. Soc. Jpn. 62, 717 (1993). 

26 L. F. Feiner, J. H. Jefferson, and R. Raimondi, Phys. Rev. Lett. 76, 4939 (1996). 

11 



T. Tohyama and S. Maekawa, Supercond. Sci. Technol. 13, R17 (2000). 

E. Pavarini, I. Dasgupta, T. Saha-Dasgupta, O. Jepsen, and O. K. Andersen, Phys. Rev. Lett. 
87, 47003 (2001). 

O. K. Andersen, A. 1. Lichtenstein, O. Jepsen, and F. Paulsen, J. Phys. Chem. Solids 56, 1573 
(1995). 

D. M. King et al, Phys. Rev. Lett. 70, 3159 (1993). 

N. P. Armitage et al, Phys. Rev. Lett. 88, 257001 (2002). 

M. Mori, T. Tohyama, and S. Maekawa, J. Phys. Soc. Jpn. 75, 034708 (2006). 

S. Ideta et al, Phys. Rev. Lett. 104, 227001 (2010). 

Y. Chen et al, Phys. Rev. Lett. 103, 036403 (2009). 

M.-H. Julien, P. Carreta, M. Horvatic, C. Berthier, Y. Berthier, P. Segransan, A. Carrington, 
and D. Colson, Phys. Rev. Lett. 76, 4238 (1996). 

H. Kotegawa, Y. Tokunaga, K. Ishida, G.-q. Zheng, Y. Kitaoka, H. Kito, A. lyo, K. Tokiwa, T. 

Watanabe, and H. Ihara, Phys. Rev. B 64, 064515 (2001). 

D. H. Kim and P. A. Lee, Ann. Phys. 272, 130 (1999). 

N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133 (1966). 



12 



